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SUMMARY 


The  application  of  the  Laplace  transformation  to  the  solution  of 
the  lateral  and  longitudinal  stability  equations  is  presented.  The 
expressions  for  the  time  history  of  the  motion  in  response  to  a 
sinusoidal  control  motion  are  derived  for  the  general  case  in  which 
the  initial  conditions,  initial  displacements  and  initial  velocities, 
are  assumed  different  from  zero.  Some  illustrative  examples  of  the 
application  of  the  Laplace  transform  to  ordinary  linear  differential 
equations  with  constant  coefficients  and  a  numerical  example  of  a 
specific  problem  are  presented  in  appendixes. 


INTRODUCTION 


Recent  developments  in  piloted  and  pilotless  aircraft,  equipped 
with  automatic  devices,  have  directed  the  attention  of  engineers  to 
the  theoretical  investigation  of  dynamic  longitudinal  and  lateral 
stability  problems  of  aircraft  designed  for  high-speed  and  high- 
altitude  flight.  In  the  past,  the  dynamic  stability  investigations 
were  usually  limited  to  the  determination  of  Routhfs  condition  for 
stability  and  for  the  calculation  of  the  roots  of  the  characteristic 
stability  equation  to  determine  the  damping  of  the  modes  of  motion 
and  the  period  of  the  oscillation.  A  more  complete  analysis  of  the 
problem  requires  the  calculation  of  a  time  history  of  the  airplane 
motion  in  response  to  a  gust  disturbance  or  in  response  to  the 
application  of  the  control  surfaces.  As  the  methods  of  classical 
analysis  (references  1  and  2)  proved  to  be  inadequate  for  this 
purpose,  new  methods  of  operational  mathematics,  representing  a 
more  powerful  tool,  were  used.  These  methods  are  known  today  as  the 
Heaviside  operational  calculus  and  the  Laplace  transformation.  The 
application  of  the  Heaviside  operational  calculus  to  the  calculation 
of  airplane  motions  is  discussed  in  references  4,  and  5-  However, 
the  Laplace  transformation  is  considered  a  more  powerful  method  than 
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the  Heaviside  operational  calculus  "because  the  initial  conditions  of 
the  problem,  initial  displacements  and  initial  velocities,  are  inher¬ 
ently  taken  into  account  by  the  Laplace  transformation,  whereas  in 
the  Heaviside  operational  calculus,  all  initial  conditions  are  zero. 

In  this  paper,  the  Laplace  transformation  is  applied  to  both  the 
longitudinal  and  lateral  stability  equations  for  the  general  case  where 
the  initial  displacements  and  initial  velocities  were  assumed  different 
from  zero.  The  operational  equations  obtained  for  this  general  case 
were  then  solved  and  the  time  history  of  the  motion  was  obtained  by  the 
Heaviside  expansion  theorem  and  by  the  inversion  theorem  for  Laplace 
transformation  -  The  Laplace  transformation  is  simple  and  effective. 

Its  principles  are  easily  understood  and  its  technique  quickly  learned. 
It  represents  a  further  development  in  operational  mathematics  because 
it  is  a  more  powerful  mathematical  tool  and  because  the  difficulties 
and  obscurities  of  the  work  of  Heaviside  are  avoided. 

A  short  historical  sketch  tracing  the  development  of  operational 
mathematics  and  its  application  to  airplane  dynamics  is  presented  in 
appendix  A. 

The  author  is  indebted  to  Mr.  Leonard  Sternfield  of  the  Langley 
Stability  Research  Division,  MCA,  for  information  and  collaboration 
he  has  contributed  in  connection  with  this  paper. 
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SYMBOLS 

chord,  feet 
span,  feet 

wing  area,  square  feet 
weight,  pounds 
mass,  slugs 
density,  slugs  per  cubic  foot 
airspeed,  feet  per  second 
time,  seconds 

nondimens ional  time  parameter  based  on  chord  (t  ^ 
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nondlmensional  time  parameter  "based  on  span  ^ 
relative  density  coefficient  based  on  chord  [ 

\pscy 

relative  density  coefficient  based  on  span 


(A 


differential  operator  with  respect  to  sc 
differential  operator  with  respect  to  s^ 

operator  in  Laplace  transformation 
root  of  stability  equation 


radius  of  gyration  about  principal  longitudinal  axis* 
feet 

radius  of  gyration  about  principal  lateral  axis,  feet 


radius  of  gyration  about  principal  vertical  axis, 
feet 

nondimens ional  radius  of  gyration  about  principal 
longitudinal  axis 

nondlmensional  radius  of  gyration  about  principal 
lateral  axis  ^ky^c^ 

nondlmensional  radius  of  gyration  about  principal 
vertical  axis  ^k£0/b^ 

angle  between  principal  longitudinal  axis  of  inertia 
and  flight  path  (fig.  l),  degrees 

angle  between  reference  axis  and  principal  longi¬ 
tudinal  axis  (fig.  l),  degrees 

angle  of  attack  (fig.  l),  degrees 
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flight-path  angle  between  path  and  horizontal 
(fig.  l),  degrees 

attitude  angle  between  reference  line  and*  horizontal 
line,  degrees  (a  +  7) 

deflection  angles  of  aileron,  elevator,  and  rudder, 
degrees 

pitching  angular  velocity,  radians  per  second  (0) 
angle  of  sideslip,  radians 
azimuth  angle,  radians 

yawing  angular  velocity,  radians  per  second  (ty) 
angle  of  hank,  radians 

rolling  angular  velocity,  radians  per  second  (fy 
increment  of  forward  velocity,  feet  per  second 

nondimens ional  increment  of  forward  velocity 

nondimens ional  radius  of  gyration  about  longitudinal 
stability  axis  (^%q^cos^tj  +  Kj^sin^riy 

nondimens ional  radius  of  gyration  about  vertical 
stability  axis  (^Kg^cos^rj  +  Kj^sin^ri^ 

nondimens  ional  product— of— inertia  parameter 


((%/  -  Kx02)Bin  n  008  n) 


rolling-velocity  parameter  (helix  angle  generated 
by  wing  tip  in  roll),  radians 


pitching-velocity  parameter,  radians 
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jawing— velocity  parameter,  radians 


dynamic  pressure 
rectangular  coordinates  (fig.  l) 
longitudinal  force,  pounds  (fig.  l) 
lateral  force,  pounds  (fig.  l) 
normal  force,  pounds  (fig.  l) 
rolling  moment,  foot-pounds  (fig.  l) 
pitching  moment,  foot-pounds  (fig.  l) 
yawing  moment,  foot-pounds  (fig.  l) 


drag  coefficient 


lift  coefficient  r^l 

V  qs  J 

longitudinal— force  coefficient 


/JL\ 

\qs) 


(: 


Y_ 

qS 


lateral— force  coefficient 

normal— force  coefficient  L 
pitching-moment  coefficient  (&) 
rolling-moment  coefficient  &) 
yawing-moment  coefficient 
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SC 7  SC 7 

c*85  =  §5; Ba  +  ^ 5r 
Cns5  =  ^ &a  + ** 


Gl,  &2,  Hi,  Hj2,  H3 


A,  A]_,  A3I 

A,  Aj,  A2,  A3J 
A,  B,  C,  D,  E 


aU  .  .  .  ^33^ 

bll  '  *  *  b33 J 

8i i (P)  •  •  •  833(H) 
h^-^(P)  ... 


functions  of  P  on  right  side  of  oper¬ 
ational  equations 


determinants 


coefficients  in  fourth-degree  characteristic 
equations 

abbreviated  coefficients  in  operational 
equation 


abbreviated  functions  of  P  in  operational 
equation 


E 


residue 


The  subscript  o  is  used  to  indicate  initial  conditions,  a  bar 
is  used  to  denote  variables  in  the  operational  equations,  and  a  dot  is 
used  to  denote  differentiation  with  respect  to  time. 


ANALYSIS 


The  purpose  of  this  paper  is  to  show  how  the  longitudinal  and 
lateral  stability  equations  can  be  solved  by  the  Laplace  transformation. 
Thus  no  attempt  is  made  to  present  a  detailed  discussion  on  the  theory 
of  Laplace  transform,  which  can  be  found  in  references  6  and  7  snd  in  the 
bibliography  presented  in  appendix  I  of  reference  6,  but  rather  to 
present  sufficient  background  of  the  theory  to  permit  a  clear  under¬ 
standing  of  its  application  to  this  particular  problem. 

If  a  function  x(t),  defined  for  all  positive  values  of  the 
variable  t,  is  multiplied  by  e-Pfc  and  integrated  with  respect  to  t 
from  zero  to  infinity,  a  new  function  x(P)  of  the  variable  P  is 
obtained;  that  is 
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*  x(P)  =  Q  ^x^dt 

This  operation  on  a  function  x(t)  is  called,  the  Laplace  transformation 
of  x(t).  The  necessary  and  sufficient  conditions  for  the  existence  of 
the  Laplace  transform  of  a  function  x(t)  are  discussed  in  reference  6. 
Let 

^  +  *1  ^TT  +  •  •  •  If  +  -  x(t)  (1) 

dtn  dt11*”1  dt 

represent  an  ordinary  linear  differential  equation  with  constant 
coefficients  a^,  .  .  •  arb-±>  an#  If  ^  substituted 

for  — ,  for  ^r-,  and  so  forth,  equation  (l)  can  be  written  in 

at  dt^ 

operational  form 

(pn  +  a^L1^  +  .  .  .  a[1__j_Dn_^  +  a^x  =  x(t)  (la) 

When  t  =  0,  the  following  initial  conditions  are  assumed: 

* 


x  =  Xq 


XI  =  3s 

1  dt 


xn— 1  = 


dn  ~*~x 
dt 


The  Laplace  transformation  of  equation  (la),  with  the  use  of  the 
letter  P  for  the  operator,  is 
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(p11  +  a]_PrL—1  +  •  .  .  an_iP  +  an)x  =  x(P) 


Term 

corresponding 

to 


+  (p^Xo  +  Pn_2x1  +  .  .  -  Pxn_2  +  xn_1) 


rPx 


+  a1(pn“2x0  +  P^x-l  +  .  .  .  PxI1_3  +  xn_2) 


Dn-1 


X 


4f 


+  an_3(p2xc  +  Pzq  +  xg) 

+  an_2(Pxo  +  xi) 

+  an— lxo 


D3x 

D2x 

Dx  (2) 


The  transform  x(P)  for  x(t)  is  taken,  from  table  I  -which  presents 
some  simple  Laplace  transforms.  A  more  complete  table  of  Laplace 
transforms  is  given  in  appendix  III  of  reference  6  and  in  appendix  A 
of  reference  7.  Appendix  B  shows  two  illustrative  examples  of  the 
application  of  Laplace  transform  to  ordinary  linear  differential 
equations  with  constant  coefficients. 


Longitudinal  Motion 

The  nondimens ional  linearized  stability  equations  for  longi¬ 
tudinal  motion  are  given  by  NACA  in  the  form: 

2HcDcu*  =  CXse&e  -f  C]CutU»  +  (cXa  +  \  %aDc)a.  +  (cx@  +  §  C x^0 

2^cDc(a  -9)  =  ChJ>e  +  C^u*  +  (c^  +  £  C^D^a  +  +  £  Cz^6 

^c^Y^c20  =  ^e+  |  CmD8eDc)6e  +  Si.**  +  (Cma  +  2  CffiDaDc)a 


>(3) 
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The  Laplace  transformation  is  demonstrated  for  the  case  in  which  the 
elevator  motion  can  be  simulated  by  the  sine  function 

$e  =  &m  Bin  asc  (4) 

where  &m  is  the  amplitude.  (In  most  cases  ^  is  assumed  to  be  1.) 
Rearranging  and  substituting  equation  (4)  into  equation  (3)  give: 

(*cDc  -  Cv)  -  (%.  +  |  Vc)“  -  ( 

=  &msin  asc 
ue 

-°vu’  +  [jS  - 1  czm)I,o  -  0Za]0  -  L( 

=  cZ5e5msin  asc 

~ Cmui u*  ~  ^1%,  +  \  ^DaP0) a  +  ^ic®Y^c^ 

■  (C»6e+IS(5eDc)6"lSlnaS'= 

In  order  to  illustrate  the  use  of  the  Laplace  transformation  for 
a  very  general  case,,  the  only  initial  condition  assumed  to  he  zero 
is  &Go  =  0;  that  is,  the  deflection  is  measured  from  its  trim 

position  before  the  maneuver  begins.  For  all  other  parameters  the 
initial  conditions  are  assumed  to  be  different  from  zero;  thus 
the  values  are  u0f,  aQ,  0Q,  and  at  sc  =  0.  The  equations  can 

then  be  written  in  general  form,  in  which  the  four  initial  disturbances 
are  combined  with  elevator  motion.  In  a  specific  problem  some  of  the 
initial  conditions  would  probably  be  zero.  For  practical  engineering 
purposes,  in  fact,  the  most  interesting  cases  are 

(l)  Disturbance  only  in  angle  of  attack  aQ  (due,  for  example,  to  a 
gust);  elevator  fixed;  all  other  disturbances  zero  (u0f  =  0Q  =  qQ  =  0) 
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(2)  Change  in  thrust ,  thus  u0f  ^  0;  elevator  fixed;  other 
disturbances  zero  (a0  =  0O  =  qQ  =  0) 

(3)  Disturbance  caused  by  elevator  motion;  other  disturbances 
zero  (no*  =  aQ  =  90  =  q0  =  0) 

Each  of  these  assumptions  greatly  simplifies  the  equations  and 
shortens  the  computations *  because  many  terms  in  equations  developed 
for  a  general  case  will  vanish. 

The  Laplace  transformation  of  equation  (3a)  can  be  written  as 
follows : 


-  %,)  »•  -  (%  +  I  V)°-  -  (%  +  5  %p)e 


=  %  +  ^c“o'  -  I  Cy-Ia0'!  -  \  %0O 


+  [(*.  -  I  CZlta)P  -  °Za]a  -  [(^0  +  |  %)”  +  %]6 

=  +  a2  +  (2|Ac  ”  2  CzI^.)a°  “  (£Pc  *  2  Czil)90 

-  (C»a  +  I  V)3  +  -  1  CmqP  -  0^5 

’  £"8e  +  i  -  S  °”Da“°  +  2^(P8°  +  «°) 


(5) 
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Equation  (5)  can  "be  expressed  in  a  shorter  form  as 

anu1  +  ai2a  +  a^e  =  Gp 
a2]_r*  +  &22a'  +  a23^  =  ^2  V 


where 


a^pu1  +  a^^ci  +  a33^  *"  ^3 


(5a) 


all 

=  2p 

cP- 

■CV 

a12 

--( 

+  I  CXDaP) 

a13 

H 

S 

+ 1  v) 

a21 

=  -C 

V 

a22 

-  (a>o  - 

■  5  %a)P  ■ 

a23 

=  — 

(2UC 

,  +  “  Cr? 

'  2  V 

a31 

=  -C 

mu« 

a  32 

-  ~(Cma 

+  1  CmDaP) 

a33 

=  2HcKY2 

■p2  1  n 

2  % 

«] 


(6) 
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Gl.=  — ^ 


a^Xg  §m  +  (2p-cUo'  — 


\  CXDaa0  -  \  CXqflo)  (Pg  +  a2) 
(p2  +  a2) 


gj(P) 

P2  +  a2 

8l(P) 

(P  +  ia) (P  -  ia) 


G-2  = 


aCz^ta  t  |”(sHo  -  |  CZDa)“°  -  (2>'c  +  J  Oz^flo  (p2  +  a2) 

W  *  a2) 


(P  +  ia) (P  —  ia) 

aCxng^5m+  t|  Ch^^^jP  +  j^Pc-^Y  (P®o  +  ^o) 

^  (p2  +  a2)" 

S3(p) 

(P  +  ia) (P  —  ia) 


2  CmDaa°  2  CJ\0' 


>](p2  +  a2) 


Now  the  system  (5)  or  (5a)  represents  three  simultaneous  algebraic 
equations  which  can  he  solved,  for  u’ ,  a,  and.  9  by  the  method,  of 
determinants .  Thus 


u* 


-  a3 


9 


(8c) 
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where  the  determinant 


A  = 


all 

a21 

a31 


a12 

a22 

a32 


al3 

a23 

a33 


(9) 


The  expansion  of  the  determinant  A  results  in  a  quart ic  equation 
in  P 


A  =  AP^  +  BP3  +  CP2  +  DP  +  E 


(10) 


which  generally  has  two  pairs  of  complex  conjugate  roots ,  namely 


g  =  -a  —  it 


1*3^1).  =  -c  ±  id 


Thus 


A  =  (P  +  a  -  ib)(P  +  a  +  ih)(P  +  c  —  id)(P  +  c  +  id)  (10a) 


The  coefficients  of  the  quartic  (equation  (10))  are 
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B  =  -  '"lu.Mc  -  V«  -  CV%=) 

+  |  (CZlnC°o  "  CmLa%)  _  2V'-  (C”4  +  C“Iu)| 

0  "  (CZU.  Cia,1  +  Cn,,CZa  _  1,CVic  +  CZJ>J.CD0  “  4C«r,U':  +  %aCxu’ 

-  c%%  -  %a%  -  Cmu>%a  "  °-a'%)  +  4(°%aCV% 

-  Cm<1CXutcZDa  +  Cm<1CZu«CXDa  ”  CXqCZu«CmDa  “  CZ(1Cmut  GXr)a 
+  CZDaCXqLCmu,)  +  2^c%2(CXu,CZa  -  CZu,C^ 

=  [|(C%%'CZq  _  %CXu»CZa  "  CmoCXu*CZr)a  +  CmDaCXu'Cze  +  %CZu»CX<x 

+  C^Zu.Cx^  -  CXaCZu,Cma  “  CX0CZU,  Cm]>x  ~  CZqCmu,CXa  “  CZ0Cmu,CXDa 
+  CZDaCmutCX0  +  Cmu.CZaCXq)  +  2^c  (Cm0CXu«  +  Cm^CZa 
+  CiruCXu,  -  C%CZ0  “  Cmu»CXa  -  Cmu.CX0)] 


cXu.(cmacZ0  “  Cffi/Za)  +  cZu«(cm0cXa  -  CmaCX0)  +  ^.(C^Xq  CZ0cXa) 
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The  other  determinants  are 

°1  a12  a13 

A1  =  G2  a22  a23  (11) 

G3  a32  a33 

all  Gl  ai3 

Ag  =  a21  G2  a23  (12) 

a31  °3  a33 


all  a12  G1 

A3  =  a2i  a22  G2  (13) 

a31  a32  G3 

When  expanded  the  determinants  can  be  written 

4  =  Glfll(P)  -  G2f12(p)  +  G3fi3(p)  (lla) 

^2  =  ^lf2l(p)  +  G2f22(p)  “  G3p23<p)  (12a) 

A3  =  G1f31(P)  -  G2f32(p)  +  G3f33(p)  (l3a) 
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where  the  minors 

fij( 

;p) 

are 

a22 

a23 

a12 

a13 

a12 

a13 

i'll  = 

f12  = 

f13  = 

a32 

a33 

a32 

a33 

a22 

a23 

a21 

a23 

ail 

a13 

all 

a13 

f2i  = 

f22  = 

f23  = 

a31 

a33 

a31 

a33 

a21 

a23 

a21 

a22 

all 

a12 

all 

a12 

f31  = 

f32  = 

f33  = 

a31 

a32 

a31 

a32 

a21 

a22 

After  expanding,  there  result 
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f12(P>  ’  {(^W*<A2)P3 

+  [A"*2  +  i(CXDa%  '  %C"iDa)]p2 

+  |(CXacmq  +  -  CX0CmDa  -  CXqCma)p 

+  (C*aS?  "  %C^)|  (!5) 

fl3(p)  ■  j[%4  %  +  “c)  +  C*l»0  -  i  CZDo)]p2 

+  |%(f  %  +  *')  +  %(2p<=  -  I  %a) 

+  §(cX]to%  -  %°Za)]p 

+  (%%  -  %°z„)4  (l6) 

— i 

f2i(P)  =  |(-2CZu^cEy2)p2 

Si* (I  % +  2^c)  _  I  %.%> 

“  (%%•  ~  Czu-S)j  (17) 
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f22<p) 


P23^P^  _ 


f3l(P) 


f32(p) 


(18) 


* 


+  ^c) 


["°v(|  %  +  2u°)  +'2D^C  +  £  %.%] 

.%)[ 


(%CV  -  cv%) 


(19) 


i  Cm_  C*  .  +  C™  .(2^  -  i  C, 


I  SA’  +  '  2  ^Da 


+  (CVSl  "  S'°Za) 


(20) 


I 


[_(CmDaMc)p2 

-  (^c0!^  “  |  CXu.CmDa  +  £  Cmu.CXDa)P 

-  fmutCXa  - 
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The  solution  of  equations  (8a),  (8h),  and  (8c),  which  will  result 
in  a  time  history  of  u' ,  "a,  and  0,  respectively,  as  a  function 
of  sc,  can  he  obtained  from  the  Heaviside  expansion  theorem  when  there 
are  simple  poles  (reference  6).  This  expansion  theorem  is  an  efficient 
method  of  finding  the  inverse  Laplace  transform  of  the  quotient  of  two 

polynomials  If,  for  example, 

F(p) 


u* 


%  =  f(p) 

A  f(p) 


where  f(p)  and  F(p)  are  polynomials  with  no  common  factors  and  the 
degree  of  f(p)  is  lower  than  that  of  F(p),  then  for  the  case  of 
simple  poles  and  distinct  roots 


where  \n  are  the  linear  and  distinct  roots  of  F(p)  set  equal  to 
zero.  The  Heaviside  expansion  theorem  is  modified  as  indicated  in 
reference  6  if  any  of  the  roots  of  F(p)  =  0  are  repeated  linear 
factors.  It  is  important  to  note  that  the  expression  for  the  Heaviside 
expansion  theorem  given  here  is  different  from  the  expression  given 
in  reference  4  because  of  the  different  transforms  of  functions  that 
are  used  in  the  Heaviside  operational  calculus  and  Laplace  transfor¬ 
mation.  However,  if  a  problem  is  consistently  followed  through  by 
either  one  of  these  two  operational  methods,  identical  solutions  will 
be  obtained. 

The  application  of  the  inversion  theorem  of  Laplace  transformation 
to  the  solution  of  equations  (8a),  (8b),  and  (8c)  by  computing  residues 
is  shown  in  appendix  C. 
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Lateral  Motion 

The  nondimens ional  linearized  MCA  standard  equations  of  motion 

are: 


Sideslip: 

2Hb(Db|3  +  V)  =  Cy&S+CY(3P+|  Oy^  +CL^+  (CL  tan  ?)*+  |  Cy^* 
Roll: 

^b(%2r)b20  +  KxZ^b2^)  =  C  l68  +  C  + 1  C  ZpDb^  + 1  C  lr V 

Yaw: 

2^b(KZ2;Db2vl,+  KZZI)t>2^)  =  Cn55  +  CnpP  +  \  Cnp1)b^  +  § 


It  is  important  to  note  that  from  the  standpoint  of  mechanics, 

Ky7.  should  be  defined  as  Kxz  =  —  (^ZQ2  —  Kx02)sin  ^  cos  ^  ln  a  right- 
hand  system  of  axes.  However,  the  definition  of  ^XZ  presented  in 

the.  symbol  list  is  used  in  the  paper  to  conform  with  recent  NACA  standard 
equations  of  motion. 

The  Laplace  transformation  is  demonstrated  for  the  case  in  which 
the  control— surface  motion  can  be  simulated  by  the  sine  function 

5  =  Sm  sin  asb  (21*) 

where  6m  is  the  amplitude.  (In  most  cases  Bja  is  assumed  to  be  1.) 
Rearranging  and  substituting  for  &  give: 
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*  (^b^b  _  CYp)p  _  +  Cl)^  +  (^b1^ 

-°2pp +  (viV  - 1  Cljh)t  ~  (|  ClT\ 

“^pP  _  (g  -  2nbKxzPb2)  0  +  ^2nb KZ2r>b2  ~  |  QopDb)  t  =  Cng&m  sin  aeb  J 

If  the  initial  conditions  (for  s^,  =  0)  are  pD,  0O,  fQ,  pQ,  and  rD  and.  the  trim  position 
is  assumed  as  zero  (thus  80  =  0),  the  Laplace  transformation  (with  P  substituted  for  p 
to  avoid  confusion  with  angular*  velocity)  for  equation  (23a)  can  he  written  as  follows:1 

(2^P  -  o^p  -  (|  0ypP  +  CL)^  +  (axbP -l^p-  CL  tan  r)t  =  0^8*  ] 

+  2>Vo  -  \  <-Tp^o  +  (2^b  -  2  Or^o 

-Cipp  +  (2ntKx2P2  -  |  CZpp)0  -  (l  CirP  -  SutKxzP2)?  =  Cigbm  -g  » 

I 

+  2^bKX2(p0o  +  Po)  “  2  CIp^o  “  2  Cir't,o  +  2^1bKXz(p|fo  +  ro) 

-CnpP  -  CrLpP  -  S^bKxzP2)^  +  - 1  Cj^p) ?  =  C^bm  -g  &  g 

-  £  QOp^o  +  2(^bKXz(P^o  +  Po)  +  2m^(p*0  +  r0)  -  |  Cnrt0 

Equations  (25)  can  he  expressed  in  shorter  form  as 

b-j  ■}  ft  +  h-^20  +  ^13^  =  ^l 
b2ip  +  b220  +  b23i  =  h2 

^31^  +  ^32^  +  ^33^  =  s3 


>  (25) 


V  (25a) 


~  2  ~  CL  tan  ?)*  =  °T58m  sin 

-  2p-bExzBb2^  =  sin  asb 


>  (23a) 


For  practical  engineering  purposes  a  simplified  case  is  of  interest  namely 
response  to  a  horizontal  gust  \|rOJ  while  pD  =  0O  =  pQ  =  Tq  -  5  =  q. 
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where 


bll  =  2bhP  “  cYp 

bi2  =  cYpP  +  CL) 

*13  =  (2^hp  -  I  cYrp  -  CL  tan  ?) 
b21  =  “CZp 


h22  =  ^bKx2P2  _  1  Cjpp) 

b23  =  C  2rP  “  2VtPxz p2) 
b3b 

b32  =  “(g  CnpP  “  2h-b%Zp2) 
b33  =  (2^Z2p2  "  \  Cnrp) 


(26) 


%  = 


cY«5ma  + 


[2hbPo  “  |  cYp0o  -  (|  cYr  ~  2ht)t0  (P2  +  a2) 

(p2  +  a2) 


hi(P) 


2 


bl(P) 

(P  +  ia) (P  —  ia) 


(27a) 
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H2  = 


V  +  [^X2^  +Pp)  cip0Q-|  C^o  +  g^Exz^Pto  +  rQ)j  (p2  +  a2) 


P2 


+  a‘ 


h2(P) 

p2  +  a2 

b2(P) 

(P  +  ia) (P  -  ia) 


(2Tb) 


H- 


_  CnsSma  +  0^0  +  2UbKlz;(P0o  +  Po)  +  ^b^fto  +  ro)  “  \  Cn^o]  (p2  +  a2) 


+  a' 


h3(p) 

P2  +  a2 

h3<P> 

(P  +  ia) (P  -  ia) 


(27c) 


Now  equations  (25)  or  equations  (25a)  represent  three  simultaneous 

algebraic  equations  which  will  he  solved  for  0,  <j,  and  f  by  the 
method  of  determinants 


j.b 

A 


0  =  ^2 
A 


A 


(28a) 


(28b) 


(28c) 
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where  the  determinants  are 


bll 

b12 

b13 

A  =  b21 

b22 

b23 

(29) 

b3i 

t32 

b33 

b12 

b13 

11 

l<f 

b22 

b23 

(30) 

H3 

b32 

b33 

bll 

H1 

bl3 

Ag  =  tgi 

H2 

b23 

(31) 

b31 

H3 

b33 

bll 

b12 

% 

A3  =  b21 

b22 

■Hg 

(32) 

b31 

b32 

H3 

If  the  values  of  equations  (26)  are  substituted  into  equation  (29)  and 
the  determinant  A  is  expanded^  a  quartic  equation  is  obtained 


A  =  AP^  +  BP3  +  CP2  +  DP  +  E  =  0 


(33) 
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which  generally  has  a  pair  of  complex  conjugate  roots  and  two  real 
roots 

■^1,2  =  -a  1-  ih  (Dutch— roll  oscillation) 

P3  =  -c  (spiral  mode) 

P4  =  -d  (rolling  subsidence) 

Thus 

A  =  (P  +  a  -  ib) (P  +  a  +  ib) (P  +  c)(P  +  d)  (33a) 

The  coefficients  of  equation  (33)  are 

A  =  8ub3(Ex2Kz2  -  Kj/) 

B  =  -2Hb2[KX2(2Kz2CY&  +  Cnr)  +  Kz2cip  -  Kxz(%  +  +  2X^0 

c  =  >1b[KX2(cYpCnr  +  ^bcnp  -  CYrcnp)  +  KZi2(cYpClp  “  CYpclp) 

-  Exz(CYpCir  +  CypCj^  +  VbCZp  -  CyrCip  -  CYpCnp) 

-  §(%%  -  %cv)| 
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D  =  |j2^bC:L  l^211  7(Kx2Cn|3  ~  KjzC  Zp)  +  Kz,2Cip  “  KXZcnpJ 

+  i[CYp(CnpC2r  ”  %Cnr)  +  %(°%%  “  CYrcnp) 

+  Cnp(cYrcZp  “  CYpCir)]  +  ^(CIpCnp  “  cnpc  Zp)J 

E  =  —■  jtan  7(0^0^  -  CjpCn^  +  C2pCnr  -  cnpC 

The  development  of  the  determinants  (30)*  ( 31 ) ,  and  (32)  gives 

Si  =  SiFn(P)  -  H2F12(P)  +  H3F13(P)  (30a) 

A2  =  — HjF2j_(p)  ^2^22^^  —  P3P 23 (  P)  (31a) 

A3  =  E±731(T)  -  H2F32(P)  +  H3F33(P)  (32a) 


where  the  minora  are 


b22 

b23 

Fi2  = 

b12 

b13 

F13  = 

b12 

b13 

Fu  = 

"l_ 

•u 

"U 

h32 

b33 

b32 

b33 

b22 

b23 

f2i  = 

b2i 

b31 

h23 

b33 

f22  = 

bll 

b31 

b13 

b33 

f23  = 

bll 

b21 

b13 

b23 

b21 

b22 

bll 

b12 

bll 

b12 

F31  = 

b31 

b32 

f32  = 

b31 

b32 

F33  = 

b21 

b22 
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The  expansions  yield 

^11  (P)  =  ^b2(%2KZ2  “  Exz2)1^  “  ^■bj%2cnr  +  KZ2CJp 

-  +  Cj|p3  +  l(ClpCnr  -  Clrcgp2  (3M 

F12(P)  =  -  CYr)  -  Ez^YpJP3 

+  [t(CYpCnr  “  CYrCnp)  -  ^(SKz^L  ~  2KXZCL  tan  7  ~  Cnp)]  p2 
+  -^(Cnp  -  %  tan  7^P  (35) 

^(P)  =  |hbKx2(-^Pt  +  Cyr)  -  MtKxzCypJp3 

+  |^(CircYp  ~  cYrCzp)  +  Pb  |pCL^Kx2tan  7  -  Exz)  +  C ip  j  P2 

+  ^(Clr  "  %tan  7)P  (36) 

P21(P)  =  -2Pb(%2Czp  -  ExzCnp)P2  +  |(CIpCnr  -  CnpCZr)p  (37) 
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F22(p)  =  j^|ab^,2:p3  _  ^b  (^^Yp  +  Cnr)T>2 

+  [§(<VYj  -  %%)  +  &*!,%]  P  -  CL  tan  7  cj-  (38) 

F23(P)  =  [vb%xzp3  +  ^b(_2EXZCYp  “  CZr)p2 

+  ~(C i  Cv  —  Cv  c,  +  4^0,  \P  -  Ct  tan  7  0,1  (39) 

2\  lr  Y0  Yr  H  b  H)  L  *0j 

F 3l(F)  =  jW^-KxzCjp  +  Kx^np)?2  +  |(cnpclp  “  cnpcZp)pJ  (^0) 

F32(P)  -  [K2kIZp3  +  Mb(-2KECYp  -  Cnp)p2  +  §(CnpCYp  -  cnf,cYp)p 

-  °L%]  ^ 

f33(p)  -  [h¥3  -  ^(aCjSoyp  +  %)p2  *  |(cYpclp  -  cYp%)p 

-  cicipj  (*•'-) 

As  indicated,  in  the  section  entitled  "Longitudinal  Motion"  the 
solution  of  equations  (28a),  (28b),  and  (28c),  which  will  result  in 
a  time  history  of  0,  0,  and  t,  respectively,  as  a  function  of  s^, 

can.  be  obtained  from  the  Heaviside  expansion  theorem  (reference  6) 
or  by  computing  residues  as  shown  in  appendix  C. 
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General  Remarks 

The  method  presented  in  this  paper  can  also  he  applied  to  cases  in 
which  the  airplane  is  equipped  with  automatic  pilots  acting  on  elevator* 
rudder*  and  ailerons.  Each  automatic  pilot  is  characterized  by  addi¬ 
tional  equations  of  motion.  Thus*  generally  speaking*  there  are  four 
simultaneous  equations  for  longitudinal  motion  and  five  simultaneous 
equations  for  lateral  motion  (two  automatic  pilots).  As  before*  the 
Laplace  transformation  is  applied  to  these  equations  and  the  problem  is 
treated  according  to  the  method  indicated  in  this  paper. 

Some  transforms  of  simulating  functions  for  control— surface  motion 
are  presented  in  appendix  D  and  simplified  methods  of  computation  of 
Laplace  transformation  are  given  in  appendix  E. 


CONCLUDING  KEMARKS 


The  application  of  the  Laplace  transformation  to  the  solution 
of  the  lateral  and  longitudinal  stability  equations  has  been  presented. 
The  expressions  for  the  time  history  of  the  motion  in  response  to  a 
sinusoidal  control  motion  were  derived  for  the  general  case  in  which 
the  initial  conditions*  initial  displacements  and  initial  velocities* 
were  assumed  different  from  zero. 


Ryan  Aeronautical  Company 

Lindbergh  Field*  San  Diego*  Calif.*  July  22*  19^9 
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APPENDIX  A 

HISTORICAL  SKETCH 


A  short  historical  sketch  on  the  development  of  the  operational 
calculus  and  its  application  to  airplane  dynamics  is  presented. 

The  fundamentals  for  the  theory  of  small  oscillations  about  a 
steady  state  of  motion  were  developed-  in  1877  by  Routh  (references  8 
and  9).  Then  as  early  as  1903  Bryan  applied  the  mathematical  equations 
of  motion  of  a  rigid  body  to  the  disturbed  motion  of  an  airplane 
(reference  10) .  In  the  following  years  the  mathematical  theory  remained 
fundamentally  in  the  form  proposed  by  Bryan,  but  the  method  of  appli¬ 
cation  was  changed  as  the  result  of  the  development  of  experimental 
research  by  the  NACA. 

During  those  years  many  scientists  were  working  on  the  problems  of 
dynamic  stability,  not  only  in  the' United  States  but  also  in  Great 
Britain,  France,  Belgium,  Germany,  and  other  countries.  In  1927  the 
equations  of  motion  were  first  expressed  in  dimensionless  form 
"by  Glauert  (reference  ll).  Jones,  Bairstow,  Zimmerman ,  and  Millikan 
(references  1,  2,  12,  13,  and  l4)  also  dealt  with  dynamic  stability 
and  their  work  is  well-known  to  the  average  engineer  in  this  country . 

The  need  for  a  means  of  describing  the  response  of  the  system 
(mathematically  similar  to  the  system  used  herein)  to  the  applied 
disturbance  was  realized  by  electrical  engineers  many  years  ago. 

In  1899  Heaviside,  impelled  by  this  need,  contributed  a  significant 
development.  In  his  electromagnetic  theory  he  originally  devised 
his  operational  calculus  for  the  solution  of  ordinary  linear  differ¬ 
ential  equations  with  constant  coefficients  and  some  of  the  partial 
differential  equations  of  applied  mathematics.  The  principles  of 
this  method  are  illustrated  in  reference  15  • 

The  significance  of  Heaviside's  contributions  were  not  recognized 
in  his  lifetime  because  of  the  inadequacy  of  the  mathematical  treatment 
and  the  obscurity  of  his  papers.  Bromwich,  making  use  of  the  theory  of 
functions  of  a  complex  variable,  explained  and  established  the  validity 
of  Heaviside’s  methods.  Bromwich’s  method  consisted  of  finding  the 
solution  of  a  given  differential  equation,  with  initial  and  boundary 
conditions,  in  the  form  of  a  complex  integral  over  a  suitable  path; 
the  choice  of  the  integrand  and  contour  is  sometimes  difficult.  Further 
research  by  Carson,  Carslaw  and  Jaeger,  and  Doetsch  (references  1 6 
to  18)  resulted  in  the  application  of  the  Laplace  transformation  to 
the  differential  equation.  Finally,  Doetsch  recognized  fully  the 
value  of  the  "inversion  theorem"  for  the  Laplace  transformation.  Thus 
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the  Laplace  transformation  is  an  important  step  forward  in  operational 
mathematics,  A  complete  treatment  of  the  subject  of  Laplace  trans¬ 
formation  can  be  found  in  references  6 ,  7*  17*  For  some  time  now 

it  has  been  recognized  that  by  applying  the  Laplace  transform  a  better 
substitute  for  Heaviside  operational  methods  can  be  obtained. 

Among  the  early  attempts  to  apply  operational  calculus  to  the 
problems  of  stability  and  control  was  a  very  fundamental  work  well- 
known  in  this  country  (reference  k) .  This  paper  deals  with  lateral 
motion  and  applies  Heaviside  operational  calculus.  Later  several 
papers  were  written  on  the  dynamic  response  of  the  aircraft  which  also 
made  use  of  the  Heaviside  method.  Some  dealt  with  tail  load  variations 
due  to  elevator  motion  (for  example,,  see  reference  19) .  Others  dealt 
with  stability  with  free  controls  (reference  20),,  stick  forces  in 
maneuvers  (reference  21 )3  and  the  behavior  of  the  airplane  equipped 
with  automatic  control  (reference  22). 
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APPENDIX  B 

ILLUSTRATIVE  EXAMPLES  OF  APPLICATION  OF  LAPLACE  TRANSFORMATION 

Example  I 

Example  I  illustrates  the  application  of  Laplace  transform  to  the 
equation 

(d2  —  3D  +  2)x  =  ea4  (t  <  0) 

The  initial  conditions  when  t  =  0  are 

X  =  Xq 
Dx  =  Xj_ 

Table  I  (transform  3)  shows  that  the  transform  of 

x(t)  =  eat 


is 


*(p) 


i 

(p  -  a) 


Thus  the  Laplace  transformation  of  the  given  equation  is  (n  =  2  when 
equation  (2)  is  applied) 


(p2  -  3p  +  2)  x  =  +  (Pxo  +  xl)  “  3x0 
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Example  II 

Example  II  applies  the  Laplace  transform  to  the  equation 
(l)3  -  2D2  +  d)x  =  4  (t  >  0) 

The  initial  conditions  when  t  =  0  are 

x0  =  1 

*1  =  2 
x2  =  -2 

Table  I  (transform  1)  gives  the  transform  for 

x(t)  =  4 


If  the  rules  of  equation  (2)  are  applied^  the  transform  can  he 
written  (n  =  3)  as 

(p3  _  2p2  +  p)5E  =  i  +  (p2x0  +  pxx  +  xg)  -  2(px0  +  Xl)  +  xQ 

=  |  +  (p2  +  2p  -  2)  -  2(p  +  2) 

+  b 
P 


+  1 
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appendix  C 


THE  INVERSION  THEOREM  OF  LAPLACE  TRANSFORMATION 


By  means  of  this  theorem  x(t)  can  he  obtained,  from  its 
transform  x(P) .  If 


then 


x( P)  =  e-Ptx(t)dt  R(P)  >  0 


x(t) 


1 

2iit 


p7+i» 

x(X)eXtdX 

U7— i« 


(Cl) 


where  7  is  a  constant  greater  than  the  real  part  of  all  singularities 
of  x(X)  and 


lim 

w  — >  oo 


The  path  (7  —  ioo.  7  +  ioo)  may  he  replaced  h j  a  circle  C  containing 
all  the  poles  of  the  integrand.  Then  x(t)  is  equal  to  2x1  times 
the  sum  of  residues  at  these  poles.  The  method  of  evaluating  the 
residues  is  shown  at  the  end  of  this  appendix. 


Example  I 


Let 


-  P3  +  P  -  k 
P^  -  2P  -  3 
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The  inversion  theorem,  with  X  substituted,  for  P  in  the  equation, 
gives 


07+i«  "3  , 

x  =  -X-  ■  +-X . ~  ^  eXtdX 

2lK  i j 7-1  co  X2  +  2X  -  3 


The  denominator  has  two  roots  -3  and  1.  The  residues  must  be  evaluated 
at  two  simple  poles  at  X  =  1  and  at  X  -  —  3  and  then  summed  in  order 
to  obtain  x. 


Example  2 


Consider  the  simultaneous  equations 


( 3p  +  2)x  +  Py  =  ^ 


which  yield 


Px  +  (4P  +  3)y  =  0 


4P  +  3 

P(P  +  1)(11P  +  6) 


(IIP  +  6)(P  +  1) 


Application  of  the  inversion  theorem  of  Laplace  transformation  gives  the 
solution 


x=_^  r+im  „  ^  ±  3  0xtd^ 

21"  U7- loo  X(x  +  1)(11X  +  6) 
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y 


1 

2in 


eXt 

(ux  +  6)(x  +  i) 


dX 


The  sums  of  residues  on  each  pole  for  x  and  y  give  the  solution  as 
a  function  of  time. 

The  inversion  theorem  for  Laplace  transformation  is  now  applied 
to  equations  (8a),  (8b),  and  (8c).  If  X  is  substituted  for  P,  new 
determinants  A,  Ax,  A2,  and  A3  are  obtained  from  A,  Aj_,  Ag,  and  A3. 

A  time  history  of  u’ ,  a,  and  9  as  a  function  of  the  parameter  sc 
is  obtained  when  the  inversion  theorem  is  applied  to  these  equations 
(the  path  angle  7  is  then  also  determined  as  6  =  a  +  7) : 


eXscdX 


(C2) 


(C3) 


(CM 


In  order  to  evaluate  the  integrals  (C2),  (C3),  and  ( Cl) ,  the  values  of 
all  residues  for  each  integral  must  be  found  and  summed.  This  procedure 
is  demonstrated  on  equation  (C2) .  The  parameter  X  is  substituted 
for  P,  equation  (10a)  is  substituted  for  A,  and  equation  (lla)  is 
substituted  for  E]_.  The  auxiliary  substitutions  for  G-^X),  G2(x), 

and  &3(X)  are  obtained  from  equations  (7a),  (7b),  and  (7c).  The 
values  for  fxl(X),  f12(X),  and  f13(X)  are  obtained  from  equa¬ 
tions  (14),  (15),  and  (l6)  and  \  is  again  substituted  for  P.  Then 
equation  (C2)  can  be  written 


Y)(ex  p  +  \)(ot  +  g  +  -^(ox:  -  q,  +  \)(t3t  -  x)(bt  +  \)  °=T->Cn  ^ 

[(Y)eej(Y)e9  +  (\)3EJ(xf3  -  (\)TTJ(X)TS]  aT+J  ~ 
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different.  (See  appendix  D.) 
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*0 

Now  u'  =  Ejl  +  E2  +  ^3  +  ^4*  Substituting  the  values  for  R-j^  R2,  R^, 
and  Rij.  yields 

•u  =  (O.OOOOI562  +  0 . 0000729^-8i)  e  0  *  739517+11  •  5635i)  sc 
+  (0.00001562  -  0 . 0000729^8i) e 0 * 739517—11 • 5635i) Bc 
+  (0.0000009546  -  11. 378551  )e^~°‘ 016^,82+0 .044l8956i)sc 
+  (O.OOOOOO9546  +  ll. 378551) e(-°*016582-0. 044189561) bc 

u'  =  e“°*73951r''Sc(o.00001562e11*5^35lBc  +  0. 00001562a-11 ’ 5635i8° 

+  0.000072948ie11‘5635l8°  -  0.000072948ie-11‘5635i80) 

+  e~° ' 01^582sc (Q  .  O000009546e“° ’ °)^18956isc 

+  0.0000009546e“^'0l^l8956iSc  -  11. 3785510°  ,0i<1|1®956iSc 

+  ll.37855ie_0*°1<J,'l8956isc) 

With  the  use  of  the  relations 

eix  +  e”1-*-  =  2  cos  x 

eix  -  e~ix  =  2i  sin  x 
the  value  of  u*  can  he  expressed  as 

u*  =e"°' 73951780 (0:00003124  cos  11.5635sc  -Q--000li)-5896  sinll.56358c) 

+  e-0, 01658280 (0.0000019092  cps  0.044l8956sc  +  22.75710  sin  0.044l8956sc) 

The  roots  must  he  computed  very  exactly  to  several  decimals;  other¬ 
wise  the  computation  hy  the  method  of  residues  does  not  check  to  zero. 


bb 
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Evaluation  of  Residues 


(a)  Simple  pole: 

The  residue  at  a  simple  pole  of  the  function  f(z)  is 


R  =  lim  (z  —  a)f(z) 
z  >a 

(b)  Multiple  pole: 

Let 

g(z)  =  (z  -  a)“f(z) 

Then  the  residue  at  the  pole  of  mth  order  is 


(C  6) 


(C7) 


Example  1  (Simple  poles).—  Let 


(C8) 


7+i« 


X3  +  X.  -  b _ 

(X2  -  2X  +  2) (x2  +  2X  -  3) 


The  denominator  has  four  roots.  There  are  four  simple  poles  (m  =  l) 
at  A>i  =  —3^  X2  =  lj,  and  =  1  ±-  i.  The  computation  of  the  residue 

at  a  simple  pole  X  =  —  3  is  illustrated.  If  the  term  under  the  integral 
sign  is  called  f  (X)dX^  then  according  to  equation  (C7) 


g(X)  =  (X  -  l)f(X) 
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Thus  the  residue  is 

E  =  (A.  -  l)(x3  +  \  _  4)ext 

I  U2  -  2X  +  2)  (x  -  l)  (x  +  3)Jx=_3 

=  (-27  -  3  -  4)e~3t 
(9  +  6  +  2) (—3  -  1) 

-  e-3t  =  i  e-3t 

-4  x  17  2 

Example  2  (Double  pole).—  Let 


(t)  =  A-  r "  + xi + 4xo)  ext( 

2ilt  Jy-ioo  (x  +  2) 2 


There  is  now  a  double  pole  (m  =  2)  at  X,  =  a  =  —2.  According  to 
equation  (C7) 

g(X)  =  (X  +  2)2f(X) 

which,  when  substituted  into  equation  (C8),  yields 


B=I1.  (X  +  2)2(Xxq  +  xq  +  4-Xp)  qU 
11  (X  +  2) 2 


[xqS^  +  t(Xxo  +  x1  +  Xxo^e^J 


X=—  2 


=  Xoe  ^  +  t(— 2Xq  +  x1  +  4x0)e  ^ 
=  [xQ  +  (2xq  +  x1)t]e“2t 
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Example  3  (Triple  pole) Let 


7+ioo 


x  - 


Jy-ico  (x2  +  a2)3 


e^dX 


There  are  triple  poles  at  X  =  ia  and  at  X  =  —la.  If  X 
then  equation  (C7)  gives 


g(X)  =  (X  -  ia)3  - 


.Xt 


(x2  +  a2) 


which,  substituted,  into  equation  (C8),  gives 


Ei-fr 


2  (X  -  ia)3XeXt 


dX  [(X  +  ia)(X  -  ia)] 


X=ia 


l6a2 


(t  ♦  jv** 


Similarly,  at  X  =  — ia 


x  =  +  I?2 


With  the  use  of  the  relations 


sin  x  = 


eix  _  Q-ix 
2i 


eix  +  Q-ix 


=  ia. 


COS  X 


2 
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then 


t 

8a3 


(sin  at  —  at  cos  at) 


X 


48 


MCA  TN  2002 


APPENDIX  D 

TRANSFORMS  OF  SIMULATING  FUNCTIONS  FOR  CONTROL-SURFACE  MOTION 

In  order  "to  include  the  control— surface  deflection  as  a  function 
of  time  (or  parameter  sc  and  s^ )  in  the  equations  of  motion,  the 

assumed  motion  must  he  simulated  (for  simplicity)  ty  some  simple  known 
f  met  ions . 

Some  examples  are  given  which  could  he  satisfactory  in  many 
practical  cases.  The  functions  are  given  here  for  five  examples  and 
their  transforms  are  given  in  table  II. 

(l)  Step  function: 


(F  =  0;  t  <  0) 
(F  =  F0;  t  >  0) 

t 


(la)  Straight  line: 


(F  =  mt  +  h) 
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(5) 


It  is  often  convenient  to  get  the  result  for  unity  of  control- 
surface  deflection  (say  one  radian) ;  then  the  result  for  any  arbitrary 
deflection  can  he  readily  computed. 
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APPENDIX  E 

SIMPLIFIED  METHODS  OF  COMPUTATION 


The  computation  of  the  time  history  of  any  parameter  can  be 
shortened  in  some  cases  if  the  form  of  the  solution  and  also  the 
value  of  the  function  and  its  first,  second,  and  third  derivatives 
at  a  given  time  are  known.-*-  It  is  preferable  to  find  these  boundary 
conditions  at  time  t  =  0  (initial  values).  On  the  assumption  that 
the  inversion  theorem  for  Laplace  transformation  gives  a  parameter  y 
in  the  form 


y  = 


l 

2  ni 


17+i 
U7— i°° 


f(X.)extdX 


(X  -  \±){\  -  \2)(\  -  x3)(x  -  A.^)X 


(El) 


this  method  will  "be  applied  to  the  longitudinal  and  lateral  motions* 

Longitudinal  motion,--  The  denominator  (of  the  stability  equation) 
has  five  roots  with  one  root  =  0.  The  other  four  roots  are 
frequently  two  conjugate  complex  pairs: 


^1  2  —  co  i.  ip 
x3,4  “  5  ±  17 

If  equation  (l)  is  assumed  to  give  the  solution  in  the  form 

y  =  CQ  +  ea^(Ci  sin  pt  +  C2  cos  pt)  +  6^(03  sin  5t  +  C4  cos  5t)  .  .  . 

(E2) 

then  CQ  =  lira  y  can  be  evaluated.  Differentiation  of  equation  (E2) 
t  — 

gives  the  derivatives  y ,  y,  and  y\ 

If  the  initial  values  of  these  functions  are  known,  that  is,  y0, 
y0,  yQ,  and  y’0,  at  t  =  0  the  constants  C-^,  C2,  C3,  and  of 

•^The  author  is  indebted  for  this  suggestion  to  Mr.  J.  M.  Debevoise. 
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equation  (E2)  can  be  obtained,  by  solving  four  simultaneous  algebraic 


equations  which  can  be  written  (for  t  =  0)  as  follows: 

yQ  =  Co  +  c2  +  C4  *  ‘  *  .  (E3) 

Yq  —  PC-^  +  cxC^  +  &c3  7^4  •  •  •  (E4) 

yQ  =  2apC1  +  (a2  -  p2)c2  +  27&C3  +  (r2  -  82)c4  .  .  .  (15) 

y'o  =  (3a2P  -  p3)Cl  +  (a3  -  3a02)c2  +  (372S  -  &3)c3 

+  (73  -  3752)c4  .  .  .  (16) 


The  value  of  *y*  can  also  be  obtained  by  plotting  y  against  t. 

y  =  eat  -  [(a2  -  P2)C1  -  2apC2]  sin  0*  +  [^aPCl  +  (a-2  -  P2)c2]  cos  Ptj* 

+  en  ^72  -  &2)C3  -  2780 ][|  sin  5t  +  [27803  +  (72  -  82)Cj[]cos  &t| 

Lateral  motion.—  For  the  case  of  lateral  motion  the  equation  has 
one  root  =  0,  one  conjugate  complex  pair,  and  two  real  roots 

^■1,2  =  a  —  iP 


X3  =  7 


X4  =  8 


The  time  history  of  any  parameter  y  can  be  written 


y  =  CQ  +  e^(C^  sin  pt  +  C2  cos  pt)  +  e^C3  +  e5^C; 
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and,  as  before,  0o  =  11m  y.  The  initial  values  of  y0,  y0,  y0, 

and  y*0  for  t  =  0  are  assumed  to  be  known.  Differentiating 

equation  (E7)  and  substituting  for  t  =  0  give  four  simultaneous 
algebraic  equations  from  which  the  constants  Cj_,  C2,  C3,  and 

of  equation  (E2)  can  be  determined*. 


yQ  =  C0  +  C2  +  C3  f  C4  .  .  .  (E8) 
yQ  =  pcx  +  aC2  +  ?c3  +  6C^  ...  (E9) 
yQ  =  2aPC1  +  (a2  -  p2)c2  +  +  52C]j.  .  .  .  (E10) 
y*o  =  (3a2P  -  P3)^  +  (a3  -  3ap2)c2  +  73C3  +  53C4  .  .  .  (Ell) 


If  desired  the  value  of  y#  can  he  obtained  by  plotting  y  against 
time 


y  =  ea^  j^(a2  —  P2)cl  —  ^aPC^jsin  pt 


+  en72 C3  +  e5*^2^ 


In  the  cases  in  which  there  are  simple  poles,  the  , method  presented 
in  reference  7  can  be  used.  This  method  is  illustrated  briefly.  It 
has  been  seen  that  any  parameter  y  of  the  equations  of  motion  (such 
as  a,  Q,  <f>,  and  so  forth)  expressed  as  a  function  of  the_ 
operator  P  can  be  written  as  a  ratio  of  two  determinants  A(P) 
and  B(P)  which  are  polynomials  in  P 


y(P) 


A(P) 

B(P) 


Thus  the  characteristic  equation  is 


B(P)  =  0 
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This  equation  is  a  polynomial  in  P,  the  highest  power  of  P  being  q; 
thus  it 'has  q  roots:  A-j_,  Ag,  •  •  •  A^.  The  inversion  theorem  gives 


y(t) 


A(t) 

B(t) 


Vs  A(Ak)  Akt 

&  B'(Xk) 


(t  £  0) 


where 


B*  =  y  B(P) 

For  multiple  poles  and  special  cases,  see  reference  7  (pp.  152 

to  169) . 
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TABLE  I 

SIMPLE  LAPLACE  TEAE5F0EMS 


Transform 

x(t) 

■NHOBMMi 

1 

1 

l 

p 

p 

t1*-!  (n  a  positive 

i 

d. 

(n  —  l) !  integer) 

pn 

3 

eat(P  >  Ee(a)) 

1 

P  -  a 

4 

sin  at 

a 

P2  +  a2 

5 

cos  at 

P 

P2  +  a2 

6 

sin  h  at(P  >  |  a| ) 

a 

P2  -  a2 

7 

cos  h  at(P  >  | a | ) 

P 

P2  —  a2 

8 

sin  at 

2a 

P 

(p2  +  a2)2 

9 

1 

. .  ^  tJl  11  Cl  0  CLO  LUb  Cl  L>  ) 

2a^ 

(p2  *  a2)2 
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TABLE  II 

TRANSFORMS  OF  SIMULATING  FUNCTIONS  FOR  CONTROL-SURFACE  MOTION 


Control— s  urf ace 
motion  determining 
x(t) 


Transform  for 
initial  condition 
80  =  0 
poo 

x(P)  =  J  e— ' ^bc(t)dt 


Transform  for 
assumption 

F0  =  1 


1 

F0(t)  step 
function 

F 

*0 

p 

1 

P 

la 

int  +  b 

(m  +  bP)/P2 

For  b  =  0,  m/P2 

2 

V** 

*0 

1 

P  4-  a 

P  +  a 

3 

F0(l  -  e^) 

F0a 

a 

P(P  +  a) 

P(P  +  a) 

4 

F0e_at(l  -  e_bt) 

F0t 

b 

(P  +  a)(P  +  a  +  b) 

(P  +  a)(P  +  a  +  b) 

5 

F0  sin  at 

F0a 

a 

P2  +  a2 

P2  +  a2 

Figure  1.-  Axes  and  notation  used.  Arrows  indicate  positive 
of  moments,  forces,  and  angular  displacements. 
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